Anti-charmed pentaquark C (3O99) from QCD sum rules 
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Abstract 

We construct QCD sum rules for the anti-charmed pentaquark © c (3099), recently reported at 
HERA. The sum rules are constructed similarly with the + (154O) sum rules using the anti- 
charmed analogue of the G + interpolating field. The strange quark and quark-gluon mixed con- 
densates, which were important in the + sum rules, are replaced by the gluon condensates whose 
contribution to the OPE is suppressed due to the heavy quark mass. Our result suggests that the 
parity of O c is positive. We identify the difference from the + sum rule, which leads to the positive 
parity in this heavy-light pentaquark system. The obtained mass is similar to the experimental 
value. 
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I. INTRODUCTION 



The exotic 6 + (1540) baryon, after its first discovery [1] and confirmations in subsequent 
experiments [2], has brought huge excitements in hadron physics. It is a narrow resonance 
containing 5 quarks, uudds, that has not been observed before. Though not confirmed, it 
is believed to be an isoscalar with spin 1/2, forming a flavor 10. NA49 Collaboration [3] 
later reported the observation of the narrow resonance, S (1862), which could be another 
member of the same multiplet 10. The observation of the anti-charmed pentaquark O°(3099) 
in the D*p invariance mass spectrum has been recently reported by HI collaboration at 
HERA [4]. It is anti-charmed analogue of the + (154O) and has the quark content uuddc. 
The pentaquark with one heavy anti-quark was first studied in Ref. [5] in a quark model. 
Then it has been studied in quark models [6] and Skyrme models [7, 8] and attracts recent 
interests [9] motivated by the current experimental activities. Thus, the pentaquarks are 
becoming a solid member in hadron spectroscopy and await a systematic compilation of 
their properties. 

One interesting model for the pentaquark is the diquark-diquark-antiquark picture of Jaffe 
and Wilczek (JW) [10]. In this picture, 6 + (1540) is composed by the constituent quarks, 
ud-ud-s. The two diquarks ud-ud are identical and the boson symmetry restricts them to 
be in a relative P-wave. Thus, the two diquarks combined with an S-wave antiquark lead 
to even parity for G + . The positive parity is also supported by the triquark-diquark picture 
[11, 12], the soliton model prediction [13], the quark potential model calculations [14] and 
quenched lattice calculation [15]. However, dynamical calculations based on QCD sum rules 
[16] or lattice calculations of Refs. [17, 18] support the negative parity of + . Therefore, 
the G + parity is an important issue to be settled and should be determined eventually 
from reaction mechanisms. Indeed, various reactions have been suggested to determine the 
+ parity. Several proposals are based on the order of magnitude of cross sections and 
polarization observables in the reactions including 'yN — > KQ + [19, 21], K + p — > n + K + n 
[20], 7n -> K~K + n [22], p + p -> E + [23] and -fN -> K*Q + [24]. 

Among various QCD sum rule calculations for + [16, 25-27], the approach proposed by 
Sugiyama, Doi and Oka (SDO) [16] is particularly interesting. Here, the interpolating field 
for + is constructed by mostly following the JW picture except for the fact that all the 
quarks are placed locally. The same current has been used in the lattice calculation [17]. 
Since QCD sum rules deal with current quarks, one can certainly form two different diquarks 
(with opposite parity) from the ud system and the boson symmetry no longer applies to 
the two-diquark system. Nevertheless, the two-diquark system still has odd parity as in 
the JW picture. The QCD sum rule of SDO has some interesting features. In the operator 
product expansion (OPE), each diquark propagates only to the diquark and the two diquarks 
with different parities do not mix each other. This feature is quite welcomed because the 
diquarks are expected to be a tightly bound system and hence they may not be easy to 
diffuse to the others in their propagations. The + properties in the SDO sum rule are 
mainly determined by the nonperturbative effects coming from the anti-strange quark. In 
particular, the negative parity, which is opposite to the JW prediction, is mainly driven 
by the quark-gluon mixed condensate, {sg s a ■ Gs), which is proportional to the average 
quark virtuality in the QCD vacuum. When the quark mass becomes heavier (like in the 
constituent quark picture), such a virtuality should become smaller and it may be possible 
to flip the parity. 

Thus, the extension to the charmed analogue C (3O99) provides an interesting test for the 
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SDO sum rule and lattice calculations [17]. Here, the charm quark is quite heavy so that the 
constituent-quark picture may fit well and the JW prediction for the parity is expected to 
be reproduced from QCD. In fact, quenched lattice calculation finds the parity of C (3O99) 
to be positive [28]. In the extension to the 6 C (3099) sum rules, there are two important 
aspects, which make this sum rule different from the SDO sum rule. First of all, since 
the charm quark is too heavy to form quark condensate, it gives non-perturbative effects 
only by radiating gluons. The quark-gluon mixed condensate {sg s cr ■ Gs), which was the 
important contribution in the + sum rule, is replaced by gluonic operators in the heavy 
quark expansion that are normally suppressed. Secondly, the charm quark mass has to 
be kept finite in the OPE, which can be done by using the momentum space expression 
for the charm-quark propagator. This is different from the light-quark sum rule where the 
calculation is performed in the coordinate space and all the quark propagators are obtained 
based on the expansion with the small quark mass. Keeping these two aspects in mind, we 
construct QCD sum rules for O c (3099) and see how they are different from the O + (1540) 
sum rule. 

This paper is organized as follows. In Section II, we introduce the interpolating field for 
the © c and show that it transforms properly under parity. Section III gives the phenomeno- 
logical side and Section IV gives the OPE side. The QCD sum rules for O c and their analysis 
are given in Section V. 



II. INTERPOLATING FIELD FOR C 

In our sum rules, we use the following interpolating field for O c , 

O c = e abk (e aef u^C l5 d f )(e b ^Cd h )TCc T k . (1) 

Here roman indices a,b, . . . are color indices, C denotes charge conjugation, T transpose. 
Note that we have introduced the 4x4 matrix T in front of the antiquark, which is to be 
determined from parity consideration below. The Cc T satisfies the charge-conjugated Dirac 
equation and represents the anti-charm quark. The diquarks, e aef uJC'j 5 df and e bgh u^Cd h , 
are isoscalar with spin 0. Due to the 75 difference, the two diquarks have opposite parities. 
To make a local interpolating field, all the quarks are defined at the same space-time. 
To determine T, we consider the parity transformation of O c , 

e> c (x') = l0 e c ( x ), x' = (t,-x). (2) 

This parity transformation must be recovered when each quark (and the antiquark) in O c 
transforms similarly, namely 

q'(x') = 7o<?(:r), (q = u,d,c). (3) 

This constraint in fact leads to the usual nucleon interpolating field [29] commonly used in 
nucleon QCD sum rules. Under this quark transformation, the two diquarks transform 

u' T (x')Cd'(x') = -u T (x)Cd(x) , 
u' T {x')C lb d'{x') = u T {x)C l5 d{x) . (4) 

The anti-quark is transformed accordingly as 

Cc' T {x') = - l0 Cc T {x) . (5) 
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Substituting these into the interpolating field Eq. (1) and demanding the consistency with 
Eq. (2), we find T = 1. This type of interpolating field with c — > s has been used to 
investigate the properties of O + (1540) [15-17]. 

III. PHENOMENOLOGICAL SIDE 

QCD sum rules for the 6 C are constructed from the following correlation function, 

11(g) =i J d A xe l<} - x (0\T(e c (x),Q c (0)\0) , (6) 

where Eq. (1) with T = 1 is used as the interpolating field. To construct the phenomeno- 
logical side, we note that the 6 C interpolating field can couple to both parities [16, 30, 31]. 
For the positive parity state, the interpolating field couples through 

(0|6 c (x)|e c (p) : P = +) = A + tf e (p)e-* x , (7) 

while for the negative parity, it couples through 

(0|e c (a;)|e c (p) :P = -) = X^ 5 U e (p)e~ ip - x . (8) 

Here, A± denotes the coupling strength between the interpolating field and the physical 
state with the specified parity. Using this, we obtain the phenomenological side of Eq. (6) 
separated into chiral even and odd parts, 

iP hcn ( g ) = -|A ± | 2 { ±me 2 + ■■■ = 4uf ca + nf en , (9) 

where the plus (minus) sign in front of me is for positive (negative) parity. The dots denote 
higher resonance contributions that should be parametrized according to QCD duality. It 
should be noted however that higher resonances with different parities contribute differently 
to the chiral-even and chiral odd parts [32]. Thus, IlP hen and H^ hen constitute separate sum 
rules. 

The spectral density is given by 

-Imn phcn (g) = ^|A±| 2 c% 2 - ml) ± m e \\±\ 2 5(q 2 - m%) + ■ ■ ■ . (10) 

7T 

We notice that the chiral-odd part has opposite sign depending on the parity while the chiral 
even part has positive-definite coefficient. Thus, the chiral-odd part from the OPE side can 
determine the parity. 

IV. OPE SIDE 

In the OPE side, we calculate the five diagrams shown in Fig. 1. To keep the charm quark 
mass finite, we use the momentum-space expression for the charm quark propagators. For 
the light quark part of the correlation function, we calculate in the coordinate-space, which 
is then Fourier-transformed to the momentum space in D-dimension. The resulting light- 
quark part is combined with the charm-quark part before it is dimensionally regularized at 
D = 4. 
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FIG. 1: Schematic OPE diagrams representing the two diquarks and anti-charm quark propagating 
from to x. The solid lines denote quark (or anti-charm quark) propagators and the dashed lines 
are for gluon. The crosses in (c) denote the quark condensate, (a) is perturbative contribution 
where the diquarks and anti-charm quark just propagate through, (b) is gluon contribution coming 
from the diquarks, (c) is the (qq) 4 contribution from the light quarks, and (d) and (e) are the gluon 
contributions from the charm quark. 



Our OPE is given by 



n ope (g) = n (a) + n (fe) + n (c) + n (d) + n (e) 



(ii) 



corresponding to each diagram in Fig. 1. The imaginary part of each diagram is calculated 

as 

-imrr^V) = -^7(qq) 4 (4 + m c ) 5( q 2 -m 2 c ) , 
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where the upper limit of the integrations is given by A = l — m 2 /q 2 . The integrations can be 
done analytically but we skip the messy analytic expressions. For the charm-quark propaga- 
tors with two and three gluons attached, we use the momentum-space expressions given in 
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Ref. [33]. In the + (154O) sum rule [16], (sg s a ■ Gs) was the important contribution. Since 
this condensate in the charm sector is replaced to (G 3 ) in the heavy quark expansion [34], we 
include (G 3 ) only from the heavy quark. The Wilson coefficients for light-quark condensates 
are found to vanish except for the (qq) A terms shown in Fig. 1(c). One can check that, for 
small m c , our OPE reproduces the corresponding OPE in Ref. [16]. The first term yields 
corresponding terms in Ref. [16], when it is truncated up to 0(m c ). The second term in 
the limit m c — > gives the same gluon condensate as in Ref. [16]. The fourth term involves 
(^-G 2 ) ^- in the limit m c — > 0. When this part is converted to the quark condensate in the 
heavy quark expansion, we find precisely the same Wilson coefficient of quark condensate 
in Ref. [16]. Note that, similarly as in the phenomenological side, the OPE has a chiral odd 
and even part, 

n ope (g)=n opc (g2)+ ^ n opc (g2) (13) 



V. QCD SUM RULES AND ANALYSIS 

QCD sum rules for O c are constructed by matching the two spectral densities in the 
Borel-weighted integral, 

/ rfgV* 2 / M2 -Im[nf c V) " ITV)] = , (i = 1, q) , (14) 

where M 2 is the Borel mass. Here, higher resonance contributions are subtracted according 
to the QCD duality assumption, which introduces the continuum threshold So- As the 
correlator contains the chiral odd and even part, we have two sum rules correspondingly, 

| A± |2 e - m |/^ 2 = / dq 2 e- q2 / M -Im[Tl° pc (q 2 )} , (15) 

Jml * 
rSo -i 

±m e |A ± | 2 e- m e/ M2 = / dq 2 e-' l2 / M -lm[U° pe (q 2 )] . (16) 

Jml 7T 

The second equation shows that the parity of O c can be determined by the sign of its 
right-hand side (RHS). We use the following QCD parameters in our sum rules [35], 

(^—G 2 ^ = (0.33 GeV) 4 , (G 3 ) = 0.045 GeV 6 , 

m c = 1.26 GeV , (qq) = -(0.23 GeV) 3 . (17) 

The value for the quark condensate corresponds to m q = 7 MeV (q = u, d) in the Gell- 
Mann-Oakes-Renner relation. The values for the quark and gluon condensates are rather 
standard in baryon or meson sum rules. For the charm quark mass, other values can be 
found in the literature, 1.1 GeV for its running mass and 1.5 GeV for the pole mass [36]. 
These fit the open charm mass in the sum rule of pseudoscalar correlator. Even a larger 
value is reported in Ref. [37]. Below, we will check the sensitivity of our result to this 
parameter as well as to the others. 

In Fig. 2, we plot the RHS of Eqs. (15) and (16) with respect to the Borel mass for various 
continuum thresholds y/S^ = 3.2, 3.4 and 3.6 GeV. As the left-hand side (LHS) of Eq. (15) is 
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M (GeV) M(GeV) 



FIG. 2: (a) The RHS of chiral even sum rule of Eq. (15) and (b) the RHS of the chiral odd sum 
rule Eq (16). The solid line is for ^So = 3.2 GeV, the dashed line for V^o = 3.4 GeV and the 
dot-dashed line for ^/So = 3.6 GeV. 

positive definite, we check whether the RHS of Eq. (15) is consistently positive. Our result in 
Fig. 2(a) indeed shows that the RHS satisfies this constraint. In Fig. 2(b), we plot the RHS 
of the chiral-odd sum rule (16). Here, we see that the RHS is positive suggesting that the 
parity of 6 C is positive. This positive parity agrees with the quenched lattice calculation [28] . 
The two observations do not change as we vary the continuum threshold from v^o — 3.2 to 
3.6 GeV. The main contribution from the OPE is from the perturbative part [Fig. 1(a)] and 
the gluon condensate coming from the light quark [Fig. 1(b)]. The contribution from (qq) 
and (G 3 ) is found to be small and therefore the uncertainties in these parameters affect the 
result marginally. Other values for m c ranging from 1.1 GeV to 1.5 GeV are found to yield 
somewhat different curves but all of them have the same sign again supporting the positive 
parity. 

Our result favoring the positive parity is consistent with the expectation from the con- 
stituent quark model picture with diquark correlations by Jaffe and Wilczek [10]. Namely, 
in the Jaffe- Wilczek picture, the two-diquark system has negative parity and, combining 
with the anticharm quark with intrinsic negative parity, the pentaquark has the positive 
parity in total. Note, in the SDO sum rule for + [16], this picture was not met by the 
large contribution from the quark-gluon mixed condensate, (sg s (r-Gs), which is proportional 
to the quark virtuality. This mixed condensate is now replaced by Fig. 1(e) through the 
heavy quark expansion, whose contribution to the OPE becomes marginal in our O c sum 
rule. This constitutes the main mechanism for yielding positive parity in this heavy-light 
pentaquark system. Therefore, we have an interesting crossover from the strange sector to 
the charm sector in the pentaquark parity. 

To check further the reliability of our sum rules, we calculate the O c mass and see if it 
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FIG. 3: The predicted + mass from chiral even sum rule (a) and chiral odd sum rule (b) deter- 
mined as explained in the text. The notations are the same as Fig. 2. 



agrees with the experimental value. The + mass is determined in two ways. 

1. We take derivative of the chiral even sum rule (15) with respect to 1/M 2 and divide 
the resulting equation by Eq. (15). As can be seen from Eq. (15), this step leads to 
—ml in the LHS. The mass is then obtained by taking negative square root of the 
resulting RHS. 

2. The same method as 1 applied to the chiral-odd sum rule, Eq. (16). 

The results are shown in Fig. 3 for the two cases. At the continuum threshold ^/S^ = 3.2 
GeV, both sum rules yield the C mass to be around 3.06 GeV and the result is practically 
independent of the Borel mass M. Changing the charm quark mass to a lower value of 1.1 
GeV leads to 3.02 GeV. For a much higher value of m c = 1.5 GeV, we have 3.14 GeV at the 
stable Borel region. So the extracted mass is slightly sensitive to the charm quark mass. The 
extracted mass is more sensitive to the continuum threshold. As we increase the continuum 
threshold, the result increases slightly also. At the large threshold of v^So — 3.6 GeV, the 
predicted mass is around 3.4 GeV. Thus, our sum rules give the O c mass that qualitatively 
agrees with the experimental value of 3.099 GeV. 

To summarize, we have constructed QCD sum rules for the recently discovered anti- 
charmed pentaquark @ c (3099). The charm quark mass is kept finite in the OPE while 
for the light quark propagators we use the coordinate space expressions obtained from the 
expansion in the small quark mass. Our sum rules suggest that the parity of O c is positive, 
which is opposite to that of O + (1540) determined from QCD sum rules. The obtained mass 
is qualitatively consistent with the experimental value of the HI Collaboration. 
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